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Abstract 

We study the renormalization group running of the tri-bimaximal mixing predicted by 
the two typical 5*4 flavor models at leading order. Although the textures of the mass matrices 
are completely different, the evolution of neutrino mass and mixing parameters is found to 
display approximately the same pattern. For both normal hierarchy and inverted hierarchy 
spectrum, the quantum corrections to both atmospheric and reactor neutrino mixing angles 
are so small that they can be neglected. The evolution of the solar mixing angle 9i2 depends 
on tan^ and neutrino mass spectrum, the deviation from its tri-bimaximal value could be 
large. Taking into account the renormalization group running effect, the neutrino spectrum 
is constrained by experimental data on 612 in addition to the self-consistency conditions of 
the models, and the inverted hierarchy spectrum is disfavored for large tan/3. The evolution 
of light-neutrino masses is approximately described by a common scaling factor. 



1 Introduction 



The neutrino physics has made great progress in the past decades. The mass square differences 
Am^Qj, Am^tj^ and the mixing angles have been measured with good accuracy [IlEllS]. A global 
fit to the current neutrino oscillation data demonstrates that the observed lepton mixing matrix 
is remarkably compatible with the tri-bimaximal (TB) mixing pattern [5], which suggests the 
following values of the mixing angles: 



The question of how to achieve TB mixing has been the subject of intense theoretical speculation. 
Recently it has been found that the flavor symmetry based on the discrete group is particularly 
suitable to reproduce this specific mixing pattern in leading order (LO). Various discrete flavor 
symmetry models have been built, please see the Refs.[.5l!B] for a review. A common feature 
of these model is to produce TB mixing at leading order, and the leading order predictions 
are always subjected to corrections due to higher dimensional operators in both the driving 
superpotential and the Yukawa superpotentials. These models provide an elegant description of 
neutrino mixing at very high energy scale, whereas the neutrino experiments are performed at 
low energy scale. In order to compare the model predictions with experimental data, one has to 
perform a renormalization group (RG) running from the high energy scale where the theory is 
defined to the electroweak scale Mz- Moreover, we note that RG effects have interesting impli- 
cations for model building, the lepton mixing angles can be magnified [7], even the bimaximal 
mixing at high energy can be compatible with low energy experiment [8]. Therefore, in a con- 
sistent flavor model building, we have to guarantee that the successful leading order predictions 
are not destroyed by the RG running corrections. The aim of this work is to analyze the RG 
corrections on the TB mixing pattern in two typical S4 flavor models [9l |10] in addition to the 
next to leading order corrections arising from high dimensional operators and to confront them 
with experimental values. We shall see that the running of the neutrino parameters is strongly 
constrained by the flavor symmetry as well, and this result holds very generally for the discrete 
flavor symmetry models. 

The 54 flavor symmetry is very interesting. From the group theory point of view [Hj , it is 
the minimal group which can produce the TB mixing in a natural way, namely without ad hoc 
assumptions. It is remarkable that we have more alternatives to realize the exact TB mixing 
than those in the A4 flavor model [91 \T0\ [T2l \T3[ [T^ I15j . In particular, the 2 dimensional 
irreducible representation of S4 group can be utilized to describe the quark sector. Moreover, 
the group 5*4 as a flavor symmetry, as has been shown for example in Refs. [HI |T71 |T8l [T9l l20]. 
can also give a successful description of the quark and lepton masses and mixing angles within 
the framework of grand unified theory (GUT). We note that S4 as a flavor symmetry has been 
investigated long ago |2H I22j. but with different aims and different results. 

The paper is organized as follows. In section [2j we briefly review the RG equations for the 
type I see-saw mechanism. Then we give a concise introduction to the Bazzocchi-Merlo-Morisi 
(BMM) 5*4 model [9] and the ^4 model of Ding [10] in section [sj where the main features of these 
models are shown. In section |4j Our results of RG effects on the neutrino mixing parameters 
for these two interesting models are presented. Finally we draw our conclusions in section [5j 
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2 Running of neutrino parameters in type I see-saw scenario 



The running of neutrino masses and lepton mixing angles is very important and has been studied 
extensively in the literature ISSllHESllMlEZlEHlEglEOlEIllSai^ in the past years. In 
particular, Antusch et al. have developed the Mathematica package REAP in Ref. |30j . which 
can solve renormalization group equations (RGE) and provide numerical values for the neutrino 
mass and mixing parameters. In this section, we present the RGEs for neutrino parameters in 
the minimal supersymmetric standard model (MSSM) extended by three singlet (right-handed) 
heavy neutrinos. The superpotential is given by 

W = D'YdQHd + U'Y^QHu + E^^YeLH^ + N'Y.LH^ + '^N'MN' (2) 

where Q and L are the left-handed quark and lepton doublets chiral superfields, respectively, C/^, 
D'^, and E'^ are right-handed up-type quark, down-type quark, heavy neutrino and charged 
lepton singlet superfields, respectively, Hu and are the well-known two Higgs doublets in 
MSSM. The Yukawa matrices Y^, Yd, Y^ and 1^ are general complex 3x3 matrices and the 3x3 
heavy neutrino mass matrix M is symmetric. Integrating out all the heavy singlet neutrinos, 
one gets the usual dimension-5 effective neutrino mass operator 

C^ = -\Kfg{Lf ■Hu){L3 -Hu) (3) 

where / and g are family indices, and the dot indicates the SU{2)l invariant contractions. After 
electroweak symmetry breaking, this operator leads to the light-neutrino masses 

mj,(/i) = sin^ /3 (4) 

where is the renormalization scale, v = 246 GeV and tan/3 = Vu/vd is the ratio of vacuum 
expectation values (VEV) of the Higgs doublets. Above the heaviest neutrino mass scale, the 
light-neutrino mass matrix reads 

= -\YJ{^,)M-\^l)Y,{ny sin2 /3 (5) 

When we evolved the energy from high energy scale down to the low experimental observation 
scale, the heavy singlet neutrinos involved in the see-saw mechanism have to be integrated out 
one by one, thus one has to consider a series of effective theories |30j . In general, the light- 
neutrino mass matrix can be written a^ 

m^ = -- / V+ 2YJm~^Yu J sin^ (3 (6) 

where the superscript (n) denotes a quantity below the nth mass threshold. In the MSSM, the 
two parts K and 2Y^ M~ Y^, evolve in the same way 

(") rp 

16^2^ = \ YlYe + YlYA x + x f yj^e + y^yA + [2 Tr(yjy, + 3y„ty„) - -gl - Qgl\ X (7) 

^We use the GUT charge normahzation for the gauge coupUng g\. 
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where t = ln(/i//io), and X stands for k or 2YJM-'^Y^. The RG equations for the Yukawa 
coupHngs Yu, Yfi, Yu, Ye and the right-handed neutrino mass matrix M are given by 
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In the fuU theory above the highest see-saw scale, the superscript (n) has to be omitted, and 
the RG equations for MSSM without singlet neutrinos can be recovered by setting the neutrino 
Yukawa couplings and the mass matrix of the singlets to be zero. Below the SUSY breaking 
scale, which is taken to be 1000 GeV in this work, we go to the standard model region. Since all 
the heavy right-handed neutrinos have already been integrated out at this scale, the neutrino 
masses are described by the effective dimension-5 operator, and the neutrino mass matrix rriy 
evolves as 
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-m,{Y}Ye) + [2Tr(3y„ty„ + W^Y^ + Y^Y^, 
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where A is the Higgs self-interaction coupling In order to calculate the RG evolution of the 
effective neutrino mass matrix, we have to solve the RG equations for all the parameters of 
the theory simultaneousljj^ At the mass threshold, we should integrate out the corresponding 
heavy neutrino and perform the tree-level matching condition for the effective coupling constant 
between the effective theories 
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3 Variants of the Two Models 

In this section, we recapitulate the main features of the 6*4 flavor model of BMM [9] and Ding 
|10) . Both models generate neutrino masses via type I see-saw mechanism, and the neutrino TB 
mixing is produced at LO. For an introduction to the group theory of 5*4 we refer to Refs. jTOlfSO] . 
the same conventions for the S4 representation matrix and Clebsch-Gordan coefficient are used 
in this work. 

^We use the convention that the Higgs self-interaction term in the Lagrangian is —^{H^H)^. 
^The running of the gauge couphngs has to be taken into account as weU, the corresponding /3 functions are 
well-known . 
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3.1 BMM ^4 model 



In this model the flavor symmetry S4 is accompanied by the cychc group and the Froggatt- 
Nielsen symmetry U{1)fn. The S4 flavor symmetry is spontaneously broken to the subgroup 
Z2 X Z2 in the neutrino sector and to nothing in the charged lepton one at leading order. This 
misalignment between the flavor symmetry breaking in the neutrino and charged lepton sectors is 
exactly the origin of the TB mixing. Furthermore, the auxiliary symmetry Z5 eliminating some 
dangerous terms, with the interplay of the continuous U{1)fn, is responsible for the hierarchy 
among the charged lepton masses. The leptonic flelds and the flavon fields of the model and 
their transformation properties under the flavor symmetry are shown in Table [T] 
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Table 1: Transformation properties of the leptonic fields and fiavons in the BMM model [S]. Note that uj is the 
fifth root of unity, i.e. u = f^'^'^l^. 



By introducing a U{l)ji symmetry, the authors in Ref. [9j have shown that the flavon flelds 
develop the following vacuum alignment at LO 
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The superpotential of the model in the lepton sector is 
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where the subscript li and I2 denote the contraction in li and I2, respectively, and dots stand 
for higher dimensional operators, which are suppressed by additional powers of the cutoff A. 
The composite X is given by 



(13) 



Taking into account the vacuum alignment in Eq. (11), the mass matrix for the charged lepton 
reads 
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where ye is the hnear combination of the ye,i contributions. The parameter u parameterizes 
the ratio v^/A, Vrj/A, v^/ A,Vip/ A and v^^/A, which should be of the same order of magnitude 
to produce the mass hierarchy among the charged fermions. The parameter t denotes the ratio 
ve/A. It has been shown that the parameters u and t belong to the range 0.01 < u,t < 0.05 [9j. 
The first term in Wu is the neutrino Dirac- Yukawa coupling, and the last two terms determine 
the mass matrix of the heavy right-handed neutrinos. Straightforwardly we have 

= ^ I 1 I yvu, MN=\b-cb + 2c -c I (15) 
^^VOIO/ \b-c -c 6 + 2c/ 

where b = 2xdVip and xt = 2xtv^. Integrating out the heavy neutrino v^, the light-neutrino 
effective mass matrix is given by the see-saw formula 





/ b+c —b+c —b+c \ 

I b{b-3c) b{b-3c) b{b-3c) \ 

-b+c -b^+4fec+3c^ b^-2&c+3c^ 

b{b-3c) b(fe2_9c2) ;,({,2_9c2) 

\ -b+c b'^-2bc+3c^ -b'^+ibc+3c^ i 

\ b{b-3c) b(fe2-9c2) 6(b2_9c2) / 



The light-neutrino mass matrix rriu can be diagonalized by 

UjniuUi, = diag{m^-^,mu2,'rnu3) (16) 
where m^i^,^2,u3 are the light-neutrino masses 



b + 3c\ 



2 2\b\ 

1 

j& + 3c| 



-3 = '^J^ (17) 



The unitary matrix Uu can be written as 

= iUrsUp (18) 

where Utb is the TB mixing matrix 

/ ^/2/3 I/VS \ 
Utb = -1/VQ l/V^ 1/V2 (19) 

V -1/Vq i/V^ -1/V2 J 

and C/p=diag(e*"^/^, e*"2/2^ g«a3/2^ matrix of phase with ai = arg(— 6 -|- 3c), a2 = arg(6) 
and as = arg(6 -|- 3c). Therefore the lepton mixing matrix is the TB mixing apart from the 
negligible corrections of order n^t^ from the charged lepton sector. We note that the right- 
handed neutrino mass matrix M^r is diagonalized by TB mixing as well, the mass eigenvalues 
are given by Mi = | — 6 -|- 3c|, M2 = 2\b\ and M3 = \b-\- 3c|. Comparing with the light-neutrino 



masses in Eq.(15), we have the interesting relation 



= ^ (20) 
2Mi ^ ' 
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The Yukawa coupling y is of 0(1), and we use |Am^j^|^/^ as the typical light-neutrino mass 
scale, then we obtain 

Mi ~ lO^^^^^GeV (21) 

The coefficients xt and Xd are expected to be of 0(1), as a consequence the VEVs and 
should be of the same order as the right-handed neutrino mass Mj. It is obvious that the model 
is rather constrained, there are only three independent parameters, which can be chosen to be 
|6| = y^v'^/ {Amy2), Z and $7. The latter two are defined according to 
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We can easily express Z and $7 in terms of the neutrino masses as 
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The above relations hold for both normal hierarchy (NH) and inverted hierarchy (IH) spectrum. 
Taking into account the experimentally measured mass difference Am^^j = ~ "^^i ^-'^d 
^'^atm — I'^vs ~ "^^i('^^2)l' have only one free parameter left, which is conveniently chosen 
to be the lightest neutrino mass. Imposing the constraint | cosOl < 1, we obtain the following 
limits for the lightest neutrino mass[^ 

m^i> 0.011 eV, NH 

m^3> 0.029 eV, IH (24) 

where the central values of Am^^j and Am^t„^ are used. We would like to stress that the 



mass squared differences are running quantities, therefore the bounds in Eq.(24) would change 
somewhat at low energy after considering the RG effects. 

The model is so predictive that we can express the Majorana phases in terms of the lightest 
neutrino mass as well. In the standard parametrization [35] , the lepton PMNS mixing matrix is 
defined by 

C12C13 S12C13 SlSe"*"^ 

UpMNS = diag(e*''% e*"^", e*"^") I -S12C23 - ci2S23Si3e"^ C12C23 - si2S23Si3e*'^ S23C13 



S12S23 - Cl2C23Sl3e*'' -C12S23 - Sl2C23'Sl3e*'^ C23C13 



X diag(e-*^l/^e-^^2/^l) (25) 

where Cij = cos 6ij, Sij = sin 6ij with 6ij £ [0,7r/2], the unphysical phases 5e, 5^ and 5r can be 
absorbed into charged lepton fields, 6 is the Dirac CP violating phase, 021 and 031 are the two 
Majorana CP violating phases, all the three CP violating phases 6, ipi and (p2 are allowed to 



vary in the range of ~ 2tt. Recalling that the leptonic mixing matrix is given by Eq.(18), we 
can identify the two CP violating phases as 

ipi = a3-ai, if2 = a3-a2 (26) 



The same parameter space is obtained by the neutrino mass sum rule method 
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As a result, we have 



-1 + 9^2 -QZsmQ 
cos ipi = , , sm (pi = , 

+ 81Z4 - 18^2 cos 20 V1 + 81Z4- 18Z2 cos2J7 

l+SZcosfi SZsinJl 

cos 030 = , sin 090 = (27) 

Vl + 9Z2 + 6Z cos Vl + 9Z2 + 6Z cos 17 

Since we can only determine cos from the neutrino mass spectrum, the Majorana phases 991, 992 
can take two sets of values corresponding to sin > and sin Q < respectively. We note that 
the Dirac CP phase is undetermined because the reactor angle is vanishing in TB mixing. The 
above successful leading order results are corrected by the NLO contributions, which consists of 
the higher dimensional operators in both the driving superpotential and Yukawa superpotentials. 
It has been shown that all the three leptonic mixing angles receive corrections of order u [9] . 



3.2 The Si model of Ding 

The total flavor symmetry of this model is 6*4 x Z3 x Z4 [lOj . It is remarkable that the realistic 
pattern of fermion masses and flavor mixing in both the lepton and quark sector have been 
reproduced in this model, and the mass hierarchies are determined by the spontaneous breaking 
of the flavor symmetry without invoking a Froggatt-Nielsen U{1) symmetry. The leptonic fields 
and the flavons of the model and their classifications under the flavor symmetry are shown in 
Table [2} where the quark fields have been omitted. 
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Table 2: The transformation rules of the leptonic fields and the flavons under the symmetry groups S4,, Z3 and 
Z4 in the Si model of Ref. [10], where lj is the third root of unity, i.e. cj = 

In this model the ^4 symmetry is broken down to Klein four and Z3 subgroups in the neutrino 
and charged lepton sector, respectively, at LO, this specific breaking scheme require flavon fields 
develop the following vacuum configuration 

{ip) = {0,V^,0), (x) = (0,F^,0) {C) = v^ 

iv) = iV^,yr,), {<P) = {V^,V^,V^), {A)=Va (28) 

We have demonstrated that this particular vacuum alignment is a natural solution to the scalar 
potential, all the VEVs (scaled by the cutoff A) K^/A, V^/A, F^/A, Vr,/A, V^/A and Va/A are 
of the same order of magnitude about 0{\'^), and this vacuum configuration is stable under the 
higher order corrections, please see Ref. jlOl for detail. Then the most general superpotential 
in the lepton sector, which is compatible with the representation assignment of Table [2| is given 
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where (..Oii, 12,2,81,32 stands for the li, I2, 2, 3i and 82 products, respectively. We note that 
and one can always set M to be real and positive by global phase transformation of the lepton 
fields, and a priori M should be of the same order as the cutoff scale A. Taking into account the 
vacuum alignment in Eq.(28), we find that the charged lepton mass matrix is diagonal at LO, 
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where ye is the result of all the different contributions of ye^- The neutrino Dirac and Majorana 
mass matrices can be straightforwardly read out as 
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where a = Vu^-x ^iid b = yu2~x- As a result, the light-neutrino mass matrix is given by 



2M 



2o? - Aab + 66^ + 2ab - 36^ + 2ab - 3b^ 
c? + 2ab - 36^ - Aab - 36^ 20^ + 2ab + 6b' 
c? + 2ab - 36^ 2o? + 2ab + 66^ c? - 4a6 - 36^ 

We can see that the mass matrix my is exactly diagonalized by the TB mixing matrix 

U^myUu = dia.g{myi,my2,iny^) 

The unitary matrix Uy is 

Uy = Utb diag(e-^"i/^e-^"2/2^g-ia3/2) 

The phases ai, 02 and 03 are closely related to the Majorana phase 

ai = arg(-(a - 3bf/M), 02 = aTg{-4a^/M), 03 = arg((a + 3bf/M) 

and the neutrino masses are given by 

myi = \{a-3bf\vl/{2M), my2 = 2\a\\l/M, m^s = \{a + 3bf\vl/{2M) 



(32) 
(33) 

(34) 

(35) 

(36) 
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It is interesting to estimate the order of magnitude for the right-handed neutrino mass M. Since 
the parameters a and b are expected to be of order A^, with this and using ^/\Am^^\ ~ 0.05 
eV as the hght-neutrino mass scale in the see-saw formula, we obtain 

M ~ 10^2/^^GeV (37) 



Similar to the analysis of section 3.1, we define 

- = i?e^* (38) 
a 

Straightforwardly we can express R and cos $ as functions of the neutrino masses 



3V m^2 



m„3 _ m„i 
m^2 rn^2 



COS ^> = "^^2 ^.2 (39) 

_|_ 2mi,3 _ -j^ 



In exactly the same way section 3.1, the Majorana phases in standard parameterization are 
determined as 

= ai - Q!3, = "2 - as (40) 

with 

(1 -9i?2)2 + 36ij2g-j^2^ ^ 12i?(l -9i?2)sin$ 



cos VI = , o^r,o„„ , Sm ipx 



1 9iZ2)2 _ 36i?2 cos2 $ ' (l-F9i?2)2_36i?2cos2$ 
-1 - 9i22cos2^> - 6i?cos$ 6i?(l 3i?cos$) sin$ ,^ , 

^°^^^= l + 9i?2 + 6iicoscI> ' ^^"^^= l + 9ii2 + 6i?cos$ ^^'^ 

Consequently, all the low energy parameters in the neutrino sector can be expressed in terms of 
the lightest neutrino mass. Imposing the condition | cos$| < 1, we get the following constraint 
on the lightest neutrino mass: 

m^x > 0.011 eV, NH 

m^3 > 0.0 eV, IH (42) 

The NLO corrections have been analyzed in detail in Ref. [10]. It is shown that both the 
neutrino masses and mixing angles receive corrections of order e ~ with respect to leading 
order result, where e parameterizes the ratio VEV/A and Ac is the Cabibbo angle. 



4 RG running effects in 6*4 flavor models 

As has been shown, the TB mixing is achieved in both BMM model and the 5*4 model of 
Ding at LO. In this section, we turn to a quantitative discussion of RG effects, and compare 
them with the NLO corrections and the experimental data. For definiteness we shall assume a 
supersymmetry breaking scale of 1 TeV, below which the SM is valid. We note that the mass 
hierarchy between top and bottom is produced via the spontaneous breaking of flavor symmetry 
in both models, and tan/3 should be small. As a result, we shall take the parameter tan/3 to 
be 10 apart except where explicitly indicated otherwise. To study the running of the neutrino 
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mixing parameters from the GUT scale to the electroweak scale, the Mathematica package 
REAP is used [30j . This package numerically solves the RG equations of the quantities relevant 
for neutrino mass and mixing, and it has been widely used for different purposes J36j. The 
package can be downloaded from http://users.physik.tu-muenchen.de/rge/REAP/index.html, 
and Mathematica version 5 or higher is required. We note that the approximate analytical 
solutions based on leading log approximation to the RG equations have been derived in Refs.pQl 
130] . which allows one to understand the generic behavior of the renormalization effects. However, 
due to enhancement /suppression factors and possible cancelations, the exact numerical solutions 
may differ considerably from those estimates. Therefore throughout this paper we adopt a 
numerical approach, exploiting the convenient REAP package. 

As has been demonstrated above, we generally need to introduce flavon fields to break the 
flavor symmetry in order to generate fermion masses and flavor mixing. In the unbroken phase of 
flavor symmetry, the flavons are active fields, therefore the corresponding RG equations should 



be modified in principle. However, the superpotentials of the models in Eq.(12) and Eq.(29) 
contain all the possible LO terms allowed by the symmetries, the invariance under the flavor 
symmetry 5*4 is maintained until we move down to the scale of the VEV of the flavon flelds, 
which is of the order of GUT scale. We conclude that the flavor structures of the models are 
preserved above the scale of the VEV of the flavon fields, the contributions of the flavon fields 
in the RG running can be absorbed by the redefinition of the model parameters [38]. In the 
following, we will discuss the RG evolution of neutrino masses and mixing parameters in both 5*4 
flavor models, stating from the initial conditions of neutrino Dirac and Majorana mass matrices 
described in section [3] at the GUT scale. In particularly, the parameter spaces are scanned. 

4.1 RG effects in the BMM models 

In this section we report results of the calculations of the RG evolution of the neutrino mixing 
parameters in the BMM model. Without loss of generality, we choose the Yukawa coupling y = 1 
for our numerical analysis. The GUT scale neutrino mass squared differences m^2 ~ ''^ui ^^'^ 
|m^3 — mf^^l are treated as random numbers in the range of 3.5 x 10"^ eV^ ~ 2.5 x lO"'^ eV^ and 
1.0 X 10-3 eV^ ~ 8.3 X 10 eV^ respectively, and the lightest neutrino mass is varied from 
the lowest bound determined by Eq.(23) or Eq.(39j) to 0.2 eV which is the future sensitivity of 



KATRIN experiment j3^ . The RG corrected neutrino mixing angles as functions of the lightest 
neutrino mass are shown in Fig. [ijfor both NH and IH spectrum]^ These plots display only the 
points corresponding to choices of the parameters reproducing Am^^j^, Arrig^j and the mixing 
angles within the 3a interval. 

We see that the lightest neutrino mass is still bounded from below, and the concrete values of 
the lower bounds are about 0.0107 eV and 0.027 eV, respectively, for the NH and IH spectrum, 
and these values are found to be almost independent of tan (3. It is remarkable that all the 
mixing parameters and Jcp are predicted to lie in a relative narrow range. For both NH and 
IH spectrum, it is obvious that the RG changes of the atmospheric and reactor angles are very 
small, the corresponding allowed regions lie within the current la bounds. In particular, the 
RG corrections to sin^ 623 and sin^ ^13 are of the same order or even smaller than the NLO 
contributions. On the other hand, the running of the solar neutrino mixing angle displays a 



^We shall show later that the mass squared difference at the GUT scale is about a factor of 1.2 ~ 3 larger than 
its low energy value in the whole spectrum. 

®The results are independent of the sign of sinf2, the reason is explained later. 
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different pattern. The RG change of sin^ ^12 is much larger than those of sin^ ^23 and sin^^is, 
which is a general property of the RG evolution |29[ [50] , consequently the deviation from its TB 
value can be large. In the case of NH spectrum and large tan /3, sin^ 612 is within the 3a limit 
only for smaller values of neutrino mass. Taking into account the lower bound on the lightest 
neutrino mass, rriui is constrained to lie in certain region, which decreases with tan f3. This point 
can be clearly seen from Fig{T] For IH spectrum, the RG effect of 612 is even larger due to the 
nearly degeneracy of mi,i and m^2- For example, for tan/3 = 10, sin^ 612 is very close or above 
the 3(T upper bound in the allowed region of m^^ , and the values of sin^ 9i2 goes completely 
beyond the 3a limit for larger tan /3. As a result, the IH spectrum is strongly disfavored for 
tan /? > 10 in the BMM model. We note that possible large deviation of solar neutrino mixing 
angle from the TB value, and small change of atmospheric and reactor angles under RG running 
are predicted as well in the Altarelli-Feruglio A4 model [38]. In Ref . |38j . the authors perform a 
general analysis of running effects on lepton mixing parameters in flavor models with type I see- 
saw, they show that, for the mass-independent mixing pattern, the running contribution from 
the neutrino Yukawa coupling Yy can be absorbed by a small shift on neutrino mass eigenvalues 
leaving mixing angles unchanged, consequently the RG change of mixing angle is due to the 
contribution coming from the charged lepton sector. This is exactly the reason why similar 
results to the A4 model are obtained here. 

The variations of Majorana phases (pi and (p2, Dirac CP violating phase S and the Jarlskog 
invariant Jqp with respect to the lightest neutrino mass are also plotted in Fig. [2} We note that 
Dirac phase 6 arises from the running effect, even though it is undetermined in the beginning. 
The initial value of Jarlskog invariant Jcp is zero due to the vanishing of the ^13 in TB mixing 
scheme, and it remains small because of the smallness of the ^13, although the value of 5 is large. 
It is remarkable that we can understand the dependence on the sign of sinJ7 exactly. At initial 
scale, the right-handed neutrino mass matrix M^r shown in Eq.(15) is complex with each other 
for sinfi > and sinJl < apart from the irrelevant overall phase, and the neutrino Yukawa 
coupling matrix can be chosen to be real. Therefore, in the case of sinO < 0, the complex 
conjugates of 1^, Yg, and k run in the same way as the corresponding quantities of sin > 
with the same initial conditions. Consequently the resulting low energy effective neutrino mass 
matrix for sin Q < is the complex conjugate of the corresponding one of sin > 0. As a result, 
the RG evolution of mixing angles and Jcp is independent of the sign of sin fi, and summation of 
the each CP phase for sinJ7 > and sinO < is equal to 27r. These results have been confirmed 
in our numerical analysis explicitly. 

Concretely the running of neutrino masses and mixing parameters with the energy scale is 
displayed in Fig. [Sjfor both NH and IH spectrum with tan/3 = 10, where the initial conditions 
for the NH and IH are chosen to be mi = 0.041 eV, Am^^i = 1.76 x lO-'^eV^, Aml^^ = 
5.85 X 10-3eV2 and ms = 0.0538 eV, Am^^j = 1.87 x lO^^eV^, Aml^^ = 5.58 x lO^^eV^ 
respectively. Reasonable values for the lower energy oscillation parameters are reached. We see 
that the deviation of the solar neutrino mixing angle 612 from the TB value can be relative large 
for the IH spectrum, the mixing angles 6*23 and ^13 and the CP phases 6, ifi and ip2 are stable 
under the RG evolution, the corresponding RG corrections are small. Since yjy^ = yH, the 
contribution from the neutrino Yukawa coupling is universal above the see-saw threshold. Then, 
only the charged lepton relevant part YeYg contributes to the change in mixing angles, and the 
evolution above the see-saw scales is essentially the same as below. This is in contrast with 
the usual situation where the neutrino Yukawa coupling plays dominant role in the running of 
neutrino mass matrix above the highest see-saw scale. Furthermore, we find that the running of 
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the neutrino mass mi,i is approximately given by a common scaling of the mass eigenvalues, this 
is the same as the situation below the see-saw scale |29t 15^ . It is remarkable that the neutrino 
mass is reduced by about 2.4 times at low energy. We note that the above results about the 
running behavior of neutrino masses and mixing parameters are very general, they almost do 
not depend on the initial conditions. 



4.2 RG effects in the 5*4 model of Ding 

It is remarkable that the heavy right-handed neutrinos are degenerate at LO, and the corrections 
to the degeneracy arising from RG running turn out to be so small that could be neglected, 
consequently the threshold effects should be very small in this case. In particular, we note 
that Y^Yiy is not proportional to the unit matrix any more, large RG effects seem possible. 



As has been demonstrated in Eq.(37), the right-handed neutrino mass M is estimated to be 
of order 10^^ ~ 10^'^GeV. Without loss of generality, we shall choose M = 10^^ GeV in the 
following numerical analysis, and we have checked that final results change very slowly with 
the parameter M. The neutrino mixing angles at electroweak scale as functions of the lightest 
neutrino mass are shown in Fig|4j it is obvious that the lightest neutrino mass for NH spectrum 
is bounded from below, and the lower bound on the lightest neutrino mass in the case of IH 
spectrum is still approximately zero. We see that the RG effects on both atmospheric and reactor 
angles are rather small, and the running of 9i2 can be large depending on tan /3 and the mass 
degeneracy. Matching 9i2 with the data already puts strong constraints on the lightest neutrino 
mass spectrum and tan /3 at the present stage, and a upper bound on the lightest neutrino mass 
is usually implied for small value of tan/3, which means that the neutrino mass spectrum can 
not be highly degenerate. In the case of tan f3 = 20, the IH spectrum is ruled out, since the value 
of sin^ 012 is much larger than its Su upper bound. While the model is within the 3a limit only 
for small neutrino mass for NH spectrum, as is displayed in Figjlj The predictions for the CP 
phases and the Jarlskog invariant are plotted in Figj5] In a similar way as section 4.1, we learn 



that the evolutions of mixing angles and Jcp do not depend on the sign of sin A, the summation 
of each CP phase for sin A > and sin A < is equal to 2tt. These points are checked by our 
detailed numerical analysis. 

The running of neutrino masses and mixing parameters with the energy scale are plotted in 
Fig|6| Being similar to the situation in the BMM model, the RG corrections to the CP phases 
6, ipi and ip2 are typically small, the corresponding curves are almost straight lines. We see that 
the neutrino mixing angles are rather stable under RG evolution except the solar angle for IH 
spectrum. The running of neutrino mass can be approximately described by a common scaling 
factor, and it reduced by about 2 times at electroweak scale. In short summary, the evolution 
of the neutrino parameters in Ding's 5*4 model is very similar to that of BMM model, although 
the textures of the mass matrices are totally different. 



5 Conclusion 

Flavor models based on discrete flavor symmetry are particularly interesting, they can produce 
the tri-bimaximal neutrino mixing (or some other mass-independent mixing patterns) at LO 
in an elegant way. It is a common feature that the LO predictions would be corrected by the 
subleading higher dimensional operators, and it have been shown that the subleading corrections 
are under control in some consistent flavor models. Since the tri-bimaximal mixing is predicted 
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at high energy scale, it is very necessary to investigate whether the RG effects would push the 
mixing parameters beyond the current allowed ranges by experimental data. 

In this paper, we have analyzed the RG running of the neutrino mass and mixing parameters 
in the BMM model and the ^4 model of Ding, both models predict tri-bimaximal neutrino 
mixing at LO, but the textures of the mass matrices are totally different. To study the running 
effects, we use the Mathematica package REAP. By detailed numerical analysis, we find that the 
evolution of neutrino mixing parameters displays approximately the same pattern in both 5*4 
models. We see that the atmospheric and reactor neutrino mixing angles are essentially stable 
under RG evolution for both NH and IH spectrum. However, the running of solar neutrino 
mixing angle depends on the neutrino mass and the parameter tan/3, and the deviation from its 
TB value could be large. After we take into account the RG effects, the neutrino mass spectrum 
is strongly constrained by the current data on 9i2, the lightest neutrino mass is bounded from 
both below and up, and the upper bound decreases with tan (3. For large tan f3 (tan /3 > 10), the 
value of sin^ 612 could be larger than its 3a upper bound for the whole spectrum in the case of 
IH spectrum. As a result, the IH neutrino mass spectrum is disfavored in the case of large tan f3. 
Moreover, we note that the running of light-neutrino masses can be approximately described by 
a common scaling factor, and they reduce by about 1.2 ~ 3 times at low energy. This effects is 
neglected in Ref.[38J. We note that the evolutions of mixing angles and Jcp don't depend on 
the sign of sinfi or sin A, and the sum of each CP phase for both sign is equal to 2tt. These 
results are confirmed both analytically and numerically. 

Finally we note that running of neutrino parameters in the Altarelli-Feruglio A4 model, 
BMM model and Ding's 5*4 model is similar to each other, although they produce tri-bimaximal 
mixing in different ways. The reason is that the neutrino Yukawa coupling only contributes 
to the running of neutrino mass, it doesn't affect the lepton mixing angles, and the change in 
mixing angles is due to the contribution from the charged lepton sector. We conclude that the 
running of mixing parameters is also severely constrained by the flavor symmetry in discrete 
flavor symmetry models. 
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Figure 1: RGE corrections to the neutrino mixing angles in the BMM S4, model with tan/3 = 10. The left column 
of the plots are the results for NH spectrum, and the right column for the IH case. In the case of sin^ 612, tan 13 = 5 
and tan /3 = 20 are considered. 
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Figure 2: RGB corrections to the CP phases and the Jarlskog invariant in the BMM S4 model with tan/3 = 10. 
The left column of the plots are the results for NH spectrum, and the right column for the IH case. 
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Figure 3: The running of the neutrino masses and mixing parameters with the energy scale in the BMM S4 model 
with tan P — 10 and Msusy = 1 TeV. The left column is the predictions for NH spectrum with mi = 0.041 eV, 
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The left column of the plots are the results for NH spectrum, and the right column for the IH case. 



21 











6*12 
013 


1 



2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
log,„(,,/GeV) 







023 
~ 




013 


Ui 



34567S9 10 
log,„(f,/OcV) 



12 13 14 15 16 



350° 
300° 
250° 
200° 
150° 
100° 
50° 
0° 



-^1 




sinA>0 

-6 


1 






_V2 





2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
log,|,(,i/GcV) 



350° 












300° 






250° 


S 






sinA>0 




200° 






150° 






100° 






50° 






0° 







2 3 4 5 6 7 8 9 10 
log,„WGcV) 



12 13 14 15 16 




350° 
300° 
250° 
200° 
150° 
100° 
50° 
0- 







sinA<0 

S 


m 


-f\ 






► — 



23456789 10 
log,„(;i/GcV) 



12 13 14 15 16 



0.12 






0.11 
0.10 
0.09 




m 


0.08 






0.07 
0.06 






0.05 


'"2 t 




0.04 


nil 




0.03 






0.02 







2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
loginWGcV) 




23456789 10 
logioWGeV) 



12 13 14 15 16 



Figure 6: The running of the neutrino masses and mixing parameters with the energy scale in the Sa model of 
Ding with tan /3 = 10 and Msusy = 1 TeV. The left column is the predictions for NH spectrum with mi — 0.0604 
eV, Am^„i = 1.43 x 10"*eV^ and Am^tm = 4.18 x 10"^eV^ The right column is for the IH case with ma = 0.0785 
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